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Abstract 

We search for lepton-flavor and lepton-number-violating r decays into a lepton (I = electron or 
muon) and two charged mesons (h, b! = -k^ or K^) using 854 fb _1 of data collected with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. We obtain 90% confidence level upper 
limits on the r — > ihh! branching fractions in the range (2.0 — 8.4) x 10~ 8 . These results improve 
upon our previously published upper limits by factors of about 1.8 on average. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) in charged lepton decays is forbidden in the Standard Model 
(SM) and highly suppressed even if neutrino mixing is taken into account. On the other 
hand, extensions of the SM, such as supersymmetry, leptoquark and many other models jl|- 
8[, predict LFV with branching fractions as large as 10~ 8 , which are accessible in current 
Z?-factory experiments. We search for neutrinoless lepton-flavor- violating t~ — > £rh + h'~ 
decays and lepton-number-violating r~ — > £ + h~h'~ decays^*', where £ is an electron or muon 
and is a charged pion or kaon. We analyse a 854 fb _1 data sample collected with 
the Belle detector [9| at the KEKB asymmetric-energy e + e~ collider [Io[ at center-of-mass 
(CM) energies at or below the T(45) and T(55) resonances. Previously, we obtained 90% 
confidence level (C.L.) upper limits on the branching fractions using 671 fb _1 of data; the 
results were in the range (3.3—16) x 10" 8 (ll| . The BaBar collaboration has also published 



90% C.L. upper limits in the range (7-48) x 10" 8 using 221 fb" 1 of data [12 . 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel- like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere joj. 

Particle identification is very important for this measurement. We use particle identi- 
fication likelihood variables based on the ratio of the energy deposited in the ECL to the 
momentum measured in the SVD and CDC, shower shape in the ECL, the particle's range 
in the KLM, hit information from the ACC, dE/dx measured in the CDC, and the particle's 
time of flight. To distinguish hadron species, we use likelihood ratios, V(i/j) = Cij (Ci + Cj), 
where Li (Cj) is the likelihood of the observed detector response for a track with flavor % 
(J). For lepton identification, we form likelihood ratios Vie) [l3| and V{p) [HJ using the 
responses of the appropriate subdetectors. 

We use Monte Carlo (MC) samples to estimate the signal efficiency and optimize the event 
selection. Signal and background event samples from generic r + r _ decays are generated by 
KKMC/TAUOLA [lj|. For signal, we generate the e + e~ — > t + t~ process, in which one r is 
forced to decay into a lepton and two charged mesons using a three-body phase space model, 
while the other r decays following SM branching ratios. Background samples from BB and 
continuum e + e~ — > qq (q = u, d, s, c) processes are generated by EvtGen [16J while Bhabha 



and two-photon processes are generated by BHLUMI [17J and AAFH [18|, respectively. In 
what follows, all kinematic variables are calculated in the laboratory frame unless otherwise 
specified. In particular, variables calculated in the e + e~ CM frame are indicated by the 
superscript "CM." 



[*] Throughout this paper, charge-conjugate modes are implied; hence, the notation t — > £hh' includes both 
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EVENT SELECTION 



We search for t + t~ events in which one r (the signal r) decays into a lepton and two 
charged mesons (h,h' = n ± or K^), and the other r (the tag r) decays into one charged 
track, any number of additional photons and neutrinos. 

For each candidate event we calculate the £hh' invariant mass (M^y ) and the difference of 
the £hh' energy from the beam energy in the CM frame (AE). In the two-dimensional distri- 
bution of Mihh' versus AE, signal events should have M^n close to the r-lepton mass (m T ) 
and AE close to zero. From MC simulations, the dominant background for the r — > fihh' 
modes is from continuum and generic r + r~ processes, while that for the r — > ehh' modes 
is from two-photon processes. Therefore, the event selection is optimized mode-by-mode 
since the backgrounds are mode dependent. This analysis includes 27% more data than 
in the previous one [ll|]. We reoptimized the selection criteria by reblinding the whole 
data set in the signal region until all selection criteria are finalized. The selection crite- 
ria are determined so that the figure of merit (FOM), defined by FOM= (Signal detection 
efficiency) / \/No. of expected background events, is maximized. Here, the expected back- 
ground is estimated using events from the simulated samples within the ±20<jM ehh , and 
±5<7a,e region on the M ihh i vs AE plane, where (TM ihh ,/AE is the resolution of M ehh i or AE, 
respectively, and details will be discussed later. 

Candidate r-pair events are required to have four tracks with zero net charge. All charged 
tracks and photons are required to be reconstructed within the fiducial volume defined by 
—0.866 < cos 9 < 0.956, where 9 is the polar angle with respect to the direction along the e + 
beam. Each charged track should have transverse momentum (p t ) greater than 0.1 GeV/c 
while each photon should have energy (E 7 ) greater than 0.1 GeV. For each charged track, 
the distance of the closest point with respect to the interaction point is required to be less 
than 0.5 cm in the transverse direction and less than 3.0 cm in the longitudinal direction. 

Using the plane perpendicular to the CM thrust axis [l9|, which is calculated from the 
observed tracks and photon candidates, we separate the particles in an event into two hemi- 
spheres. These are referred to as the signal and tag sides. The tag side is required to include 
a single charged track while the signal side is required to contain three charged tracks. We 
require one of three charged tracks on the signal side to be identified as a lepton. The 
electron and muon identification criteria are V(e) > 0.9 with momentum p > 0.6 GeV/c and 
V(fi) > 0.95 with p > 1.0 GeV/c, respectively. In order to take into account the emission 
of bremsstrahlung photons from the electron, the momentum of the electron candidate is 
reconstructed by adding to it the momentum of every photon within 0.05 radians of the 
electron track direction. The electron (muon) identification efficiency is 89% (81%) while 
the probability to misidentify a pion as an electron (a muon) is 0.7% (1.2%). 

Charged kaons are identified by a condition V(K/ir) > (0.6 — 0.9) for each mode, as 
shown in Table [J while charged pions are identified by the requirement V(K/tt) < 0.6. 
Furthermore, we apply a proton veto for kaon candidates, V(p/K) < 0.6, to reduce protons 
incorrectly identified as kaons. The kaon (pion) identification efficiency is around 80% (88%) 
while the probability to misidentify a pion (kaon) as a kaon (pion) is below 10% (12%). In 
order to suppress background from photon conversions (i.e. 7 — > e + e~), we require each h 
or h' candidate track to have V(e) < 0.1. Furthermore, we apply the condition V(fi) < 0.1 
to suppress two-photon background from e + e~ — > e + e~/i + /i~. 

Because a signal decay has no neutrino, the missing momentum is entirely due to neutrinos 
emitted from the tag side. The missing momentum p m i SS is defined as a difference between 
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TABLE I: Selection criteria for kaon identification V(K/ir) and magnitude of thrust (T). 



Mode 


V(K/ir) 


T 


T 


— > fiirir 




0.90 < T < 0.98 


T 


— > fiKir 


> 0.9 


0.92 < T < 0.98 


T 


fiKK 


> 0.8 


0.92 < T < 0.98 


T 


— > eirir 




0.90 < T < 0.97 


T 


— > eKir 


> 0.8 


0.90 < T < 0.97 


T 


eKK 


> 0.6 


0.90 < T < 0.98 



the sum of the e + and e~ beam momenta and the vector sum of the momenta of all tracks 
and photons, where photons are measured in the ECL as clusters to which no charged 
tracks are associated. Since the direction of p m i SS should lie within the tag side of the 
event, the cosine of the opening angle between p m ; ss and the charged track on the tag side 
in the CM system, cos#° a g_ miss , is required to be greater than zero. If the track on the 
tag side is a hadron, we also require cos0°^_ miss < 0.85 for the r — > fihh' modes. This 
requirement reduces continuum background with missing energy due to neutrons or K^s 
since the masses of the neutron and K\ are larger than that of the neutrino. We also require 
that cos0°^_ miss < 0.96 for r — > ehh' modes. This requirement reduces Bhabha, inelastic 
vector meson-photoproduction, and two-photon background, since these processes produce 
electron in the tag-side in many cases and the electron can produce the radiated photons 
which result in missing momentum if they overlap with the ECL clusters associated with the 
tag-side track j20[. In addition, in order to ensure that the missing particles are neutrinos 
rather than photons or charged particles that leave the detector acceptance, we impose 
requirements on p m i SS : We require that |p ^ i ss |, the magnitude of the transverse component 
of Pmiss, be greater than 0.5 GeV/c (0.7 GeV/c) for the r — > fihh 1 (ehh') modes, and that its 
direction point into the fiducial volume of the detector. For the r _ — > e~7r + 7r~ mode only, 
we apply the tighter selection requirement Ip^issl > 1-5 GeV/c. 

To reject continuum, Bhabha and fi + fi~ background, we require the magnitude of the 
thrust (T) to be in the ranges given in Table [B 

To suppress the BB and continuum background, we require that the number of photons 
on the tag side be n^ AG < 2 and n^ AG < 1 for decays with hadronic and leptonic tags, 
respectively. A leptonic tag is defined as V(e) > 0.1 or V(fi) > 0.1 while a tag is hadronic if 
the leptonic requirements are not satisfied. We allow at most one additional photon on the 
signal side. The reconstructed mass of the tag side, combining the photons with the charged 
track (assumed to be a pion mass) from the tag side, m tag , is required to be less than 1.0 
GeV/c 2 to reduce the continuum background. 

Photon conversions can result in large backgrounds when an electron is reconstructed as 
a hadron and still passes the electron veto. For the r — > ehh' modes, the e~h + and h'~h + 
invariant masses for the r~ — > e~h + h'~ modes (e + h~ and e + h'~ for the t~ — > e + h~h'~ 
modes), assigning the electron mass to both tracks, are required to be greater than 0.2 
GeV/c 2 to reduce photon conversion and other backgrounds. 

For the r — > fihh' modes, a real muon track can come from a kaon decaying in the CDC 
(K ± — > f^Vp)- Therefore, we apply a kaon veto, V(K/ir) < 0.6, for muon candidate tracks 
if the tag side track is a hadron (see Fig. [1] (a)). Another significant continuum background 
is from di-baryon production with a proton on the tag side. To suppress this background, 
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FIG. 1: 



(a) V(K/ir) for muon tracks, (b) V(p/ir) and (c) V(p/K) for hadronic tags, for r~ — > 
fi~ir + 7r~ candidate events. Signal MC (r _ — > /i _ 7r + 7r _ ) distributions are normalized arbitrarily 
while the background MC distributions are normalized to the data luminosity. The selected regions 
are indicated by arrows. 



we apply a proton veto, V(p/n) < 0.6 and V{p/K) < 0.6, as shown in Figs. [T] (b) and (c). 

For the r _ — > fi^7i + ii~ mode, we reduce the e + e~ — > /i + yU~7 background (with the 
photon converting into an electron-positron pair) by requiring that the invariant mass of the 
pion pair, evaluated by assuming the electron mass for both tracks, exceed 0.2 GeV/c 2 . In 
addition, we require the momentum of a muon in the CM system be less than 4.0 GeV/c 
and cos ^° a g_ m i ss < 0.97 if the track is a muon candidate with V(fi) > 0.1. 

For the r — > £ttK modes, the background in the signal region is from r~ — > 7t~7t + 7t~v t 
events in which both the kaon and lepton candidates are misidentified. To reduce this 
background, we require the invariant mass of the three charged tracks on the signal side, 
Mn-n-n, reconstructed by assigning the pion mass to the tracks, be larger than 1.52 GeV/c 2 
(see FigHD. 

Finally, to suppress backgrounds from generic r + r~ and continuum processes, we apply 
a selection based on the magnitude of the missing mass squared m^^. The variable is 
defined as E^ iss — p"^^, where E miss = E totll i — E vis , E tota i is the sum of the beam energies and 
E vis is the total visible energy. We apply different selection criteria depending on the type 
of a one-prong tag: there should be two emitted neutrinos (a single missing neutrino) if the 
tagging track is leptonic (hadronic). For the r —> ehh', finn and \iKK modes, we impose 
the requirements —1.5 (GeV/c 2 ) 2 < m 2 ^ < 1.5 (GeV/c 2 ) 2 for the hadronic tag and —1.0 
(GeV/c 2 ) 2 < m^ss < 2.5 (GeV/c 2 ) 2 for the leptonic tag. For the r —> jiiiK modes, where the 
residual background after all selections is larger than in other modes, we require the following 
relation between p m i SS and m^^: Pmiss > —8.0 x m^ ss — 0.5 and p m i SS > 8.0 x m^ ss — 0.5 
for the hadronic tag and p miss > —9.0 x m 2 ^ + 0.4 and p m i SS > 1-8 x m 2 ^ — 0.4 for the 
leptonic tag; here p m i SS is in GeV/c and m miss is in GeV/c 2 (see Fig. [3]). Typically, 75% 
of the generic r + r _ background is removed by these m^ iss requirements while 75% of the 
signal events are retained. 



SIGNAL AND BACKGROUND ESTIMATION 

For all modes, the M^y and AE resolutions are obtained from fits to the signal MC 
distributions, using an asymmetric Gaussian function that takes into account initial-state 
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FIG. 2: Invariant mass distribution of three charged tracks on the signal side with the pion mass 
assigned to each track (M niTn ) for t~ — > fi~n + K~ candidate events. Signal MC (r _ — > fj>~ir + K~) 
distributions are normalized arbitrarily while the background MC distributions are normalized to 
the data luminosity. The selected regions are indicated by the arrow. 



radiation. These Gaussians have widths as shown in Table [TT1 

To evaluate the branching fractions, we use elliptical signal regions that contain ±3o" of 
the MC signal satisfying all selection criteria. The ellipse is denned as 

((M thh , - M,V) cos6 - (AE - AE°) sing) 2 | ((M thh , - M° hh ,) sing + (AE - AE°) cosfl) 2 
(3v Mehh ,) 2 (3a AE ) 2 

(1) 

where M® hh , and AE° are the coordinates of the center of the ellipse and &M lhh , (&ae) is 
the average of ; and offi , ( a A^E ano ^ ^ae) in Table [TO These elliptical regions are 
determined by scanning M® hh ,, AE° and 9 to maximize the significance in MC simulation 
and obtain the highest sensitivity. Table [TTT1 summarizes the signal efficiencies for each mode. 
We blind the data in the signal region until all selection criteria are finalized to avoid bias. 

Figures H] and show scatter plots for data and signal MC samples distributed over a ±20cr 
rectangular region in the M^y — plane for the r — > \ihh! and ehh' modes, respectively. 
For the r — > fiirir modes, the dominant background is from continuum processes while small 
background contributions come from generic r + r~ events in the AE < GeV and M^hh' < 
m T region, which are combinations of a misreconstructed muon and two pions. For the 
r — y fiKK modes, the dominant background sources are continuum and r + r~ processes 
with a pion misidentified as a kaon. For the r — > finK modes, the dominant background is 
from generic r + r~ decays with combinations of a misidentified muon, a misidentified kaon 
and a real pion from r~ — > 7r~7r + 7r _ z/. If a pion is misidentified as a kaon, the reconstructed 
mass from generic r + r~ background can be greater than the r lepton mass since the kaon 
mass is greater than that of the pion. For the r — > ehh' modes, the dominant background 
originates from two-photon processes, while the background from continuum and generic 
t + t - processes is small due to the low electron fake rate. 

In order to estimate the background in the signal region, we use the number of data events 
observed in the sideband region inside a ±5<7a£ band around AE = GeV but excluding the 



S 





FIG. 3: Scatter plots of p m i ss vs. m^ aiss for r~ — > fj,~7r + K~ modes, (a) and (c) show the signal 
MC (r~ —> fj>~ir + K~) and generic t + t~ MC distributions, respectively, for the hadronic tags while 
(b) and (d) show the same distributions for the leptonic tags. Selected regions are indicated by 
lines. 



signal region. The events in this band are projected onto the Mn^y axis. For the r — y utttt 
modes, the extrapolation to the signal region is performed by fitting the sideband of the 
projected M^ Tn distribution using the sum of an exponential and a first-order polynomial 
function for generic rr and continuum, respectively. For the r — > ehh', fiKK and finK 
modes, the background remaining after all selections is small, and extrapolation to the 
signal region assumes that the background is linear as a function of M^. The systematic 
uncertainty due to the estimation of the expected background includes the contributions due 
to the statistics of the background sample and the shape of the background distribution. 
By varying the assumptions about the background shape, we verify that this effect on the 
systematic uncertainty is less than 20% and is smaller than the background statistical error. 
The signal efficiency and the number of expected background events with its uncertainty, 
obtained by adding statistical and systematic uncertainties in quadrature for each mode, 
are summarized in Table IIII1 After estimating the background, we open the blinded regions. 
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TABLE II: Summary of Mghh' an d AE resolutions (<r 



high/low 



(MeV/c 2 ) and aJ| h/low (MeV)). Here 



0-high (u 10 ") is the standard deviation on the higher (lower) side of the peak. 



T 

T~ 
T~ 
T~ 
T~ 
T~ 
T~ 

T~ 
T~ 
T~ 
T~ 
T~ 



Mode 



> [A TT^TT 

> /i + 7T~7T~ 

> e _ 7r + 7r _ 

> e + 7r~7r~ 
fj,-K + K~ 
H+K~K~ 
e~K + K~ 

■ IJL~1T + K~ 

> e~ir + K~ 

■ n~ K + 7T~ 

> e~K + Ti~ 

■ n + K~ir~ 



a 



high 
M fhh' 


low 


high 
°AE 


low 
a AE 


5.3 


5.8 


14.1 


20.1 


5.4 


5.7 


14.2 


20.1 


5.7 


6.2 


14.3 


22.0 


5.6 


6.3 


14.4 


22.3 


3.4 


3.6 


12.9 


17.2 


3.4 


3.3 


12.9 


17.3 


4.4 


4.4 


13.3 


19.8 


3.8 


4.2 


12.4 


19.9 


4.4 


4.8 


14.2 


18.8 


4.8 


5.5 


14.0 


21.0 


4.6 


5.1 


14.3 


18.7 


4.9 


5.4 


13.9 


21.2 


4.5 


4.7 


14.7 


18.6 


5.0 


5.4 


14.0 


21.2 



We observe one candidate event for each of the r" —> /i + 7r~7r~ and /i~7r + i\"~modes, and no 
candidate events for the other modes. The numbers of events observed in the signal region 
are consistent with the expected background levels. 

The dominant systematic uncertainties for this analysis come from the resolution in M^h' 
and AE and from particle identification. We estimate the uncertainties from resolutions of 
Mihh' and AE due to the difference between data and MC samples to be (3.7 — 4.8)%. The 
uncertainties due to lepton identification are 2.2% and 1.9% from modes with an electron and 
a muon, respectively. The uncertainties due to hadron identification are 1.3% and 1.8% for 
pion and kaon candidates, respectively. The uncertainty due to the charged track finding is 
estimated to be 0.35% per charged track; therefore, the total uncertainty due to the charged 
track finding is 1.4%. The uncertainty due to integrated luminosity is estimated to be 1.4%. 
The uncertainties due to the trigger efficiency and MC statistics are negligible compared 
with the other uncertainties. All these uncertainties are added in quadrature, and the total 
systematic uncertainties for all modes are (5.5—6.7)%. 



UPPER LIMITS ON THE BRANCHING FRACTIONS 



Since no statistically significant excess of data over the expected background in the signal 
region is observed, we set upper limits on the branching fractions of the r — > £hh' modes 



using the Feldman- Cousins method [2lJ. The 90% C.L. upper limit on the number of 
signal events including systematic uncertainty (sgo) is obtained using the POLE program 
without conditioning 22| based on the number of expected background events, the number 
of observed events and the systematic uncertainty. The upper limit on the branching fraction 
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TABLE III: Summary of upper limits for each mode. The table shows the signal efficiency (e), 
the number of expected background events (N-qg) estimated from the sideband data, the total 
systematic uncertainty (cr sys t), the number of observed events in the signal region (A" b s ), 90% C.L. 
upper limit on the number of signal events including systematic uncertainties (sgo) and 90% C.L. 
upper limit on the branching fraction (B) for each individual mode. 



T 

T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 
T~ 



Mode 



> e _ 7r + 7r _ 

e~K+K- 

■ n~n + K~ 

> e~ir + K~ 

■ n~K + ir~ 

> e~K + ir~ 



a. 



syst 



N ohs s 90 B (10- 



5.83 0.63 ±0.23 5.7 1.87 2.1 

6.55 0.33 ±0.16 5.6 1 4.01 3.9 
5.45 0.55 ±0.23 5.7 1.94 2.3 

6.56 0.37 ±0.19 5.5 2.10 2.0 
2.85 0.51 ±0.19 6.1 1.97 4.4 
2.98 0.25 ±0.13 6.2 2.21 4.7 
4.29 0.17 ±0.10 6.7 2.29 3.4 
4.64 0.06 ±0.06 6.5 2.39 3.3 
2.72 0.72 ±0.28 6.2 1 3.65 8.6 
3.97 0.18 ±0.13 6.4 2.27 3.7 
2.62 0.64 ±0.23 5.7 1.86 4.5 
4.07 0.55 ±0.31 6.2 1.97 3.1 
2.55 0.56 ±0.21 6.1 1.93 4.8 
4.00 0.46 ±0.21 6.2 2.03 3.2 



(B) is then given by 



B(t -> ihti) < -4? 2 -, (2) 



where N TT is the number of t + t pairs, and e is the signal efficiency. The value N TT = 782 x 
10 6 is obtained from the integrated luminosity and the cross section of r-pair production, 



which is calculated in KKMC [2j to be a TT = 0.919 ± 0.003 nb and a TT = 0.875 ± 0.003 nb 
for 782 fb" 1 at T(4S') and 72 fb -1 at T(55'), respectively. Table HLT1 summarizes the results 
for all modes. The upper limits for the r — > ehh' modes are in the range (2.0 — 3.7) x 10 -8 
while those for the r — > fihh' modes are in the range (2.1 — 8.6) x 10~ 8 . These results improve 
upon our previously published upper limits ll[ by factors of about 1.8 on average. This 



improvement results both from using a larger data sample and from the introduction of an 
improved rejection of specific backgrounds, such as di-baryon production in the continuum 
for the r — > fihh' modes, and improved kinematic event selections, for example, the M Tn7T 
requirement for the r — > inK modes. 



SUMMARY 

We have searched for lepton-flavor and lepton-number-violating r decays into a lepton 
and two charged mesons using 854 fb _1 of data. We obtain 90% C.L. upper limits on the 
branching fractions of r — > ehh' in the range (2.0 — 3.7) x 10 -8 and upper limits on r — > fihh' 
in the range (2.1 — 8.6) x 10 -8 . These results improve upon our previously published upper 



11 



a 02 

LU 

< 



-0.4 



(a) x"— >u"7t + 7t~ 



0.2 > /* t 



1.7 



5 0.2 



-0.4 



(b) x^]x + :i 

• • 

• ... 


in 

• 
• 

- ■ • 






v ■ 


- ■ * 


->.»i 
• i. 




• 






1.7 



M lnm (GeV/c' 



> 

<" 2 



-0.2 



(c) x"^u~K + K~ 



I : L_ 

1 .75 1 .8 



M^CGeV/cl 



§ 0.2 - (d) x"->n + K"K" 



-0.2 - 



1.75 1. 



M^JGeV/c' 



> 

cu 

13 0.2 



- 



-0.2 - 



(e) X ^>U~7t + K~ 



°7« 



_i_ 
1.7 



M^ K (GeV/c^ 



s °- 2 

UJ 

< 



(f) T — >|J. K + 7I 



o - 



-0.2 - 



1.7 



1.8 

M^ K (GeV/c 2 ) 



> 

13 0.2 



-0.2 



(g) X— >n 7t"K" 



1.7 



M^ K (GeV/c 2 



FIG. 4: Scatter plots in the Mehh' ~~ plane within a ±20cr region for the (a) r~ — > ^~7r + 7r~, (b) 
r~ — )• /i + 7r~7r~, (c) r~ — )• fj,~K + K~ , (d) t~ — >■ [i + K~ K~ , (e) r~ — > fi~Tr + K~ , (f) r~ — > K + ir~ , 
and (g) t~ —7- /i + 7r~iv~ modes. The data are indicated by the solid circles. The filled boxes show 
the MC signal distribution with arbitrary normalization. The elliptical signal regions shown by 
solid curves are used for evaluating the signal yield. The region between two horizontal solid lines 
excluding the signal region is used as a sideband. 



limits by factors of about 1.8 on average. These more stringent upper limits can be used to 
constrain the space of parameters in various models of new physics. 
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